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Snake states are trajectories of charge carriers curving back and forth along 
an interface. There are two types of snake states, formed by either inverting 
the magnetic field direction or the charge carrier type at an interface. Whereas 
the former has been demonstrated in GaAs-AlGaAs heterostructures, the latter 
has become conceivable only with the advance of ballistic graphene where a 
gapless p-n interface governed by Klein tunneling can be formed. Such snake 
states were hidden in previous experiments due to limited sample quality. Here 
we report on magneto-conductance oscillations due to snake states in a ballistic 
suspended graphene p-n-junction which occur already at a very small magnetic 
field of 20 mT. The visibility of 30% is enabled by Klein collimation. Our finding 
is firmly supported by quantum transport simulations. We demonstrate the high 
tunability of the device and operate it in different magnetic field regimes. 

A magnetic field fundamentally modifies the transport properties of an electronic con¬ 
ductor by acting on its charge carriers via the Lorenz force. The most prominent magne¬ 
totransport effect is the quantum Hall effect in a two-dimensional electron gas (2DEG). A 
strong perpendicular magnetic field forces the charge carriers into one-dimensional conduc¬ 
tion channels that flow along the edges of a sample. At moderate magnetic fields however 
electron trajectories can be understood in a quasiclassical picture where the Lorenz force 
bends charge carriers into cyclotron orbits. At the boundary of a conductor charge car¬ 
riers are not localized and can propagate via so-called skipping orbits. Magnetic focusing 
experiments [1, 2] represent a direct proof of the skipping orbit picture. In such experiments 
an increase of conductance is observed if the distance between two contacts is an integer 
multiple of the diameter of a cyclotron orbit. One condition for the observation of such 
trajectories is ballistic transport over the relevant device dimensions. This has limited the 
observation of skipping orbits to the cleanest available semiconductor samples. 

Since 2004 graphene as a new two-dimensional conductor has moved into the focus of con¬ 
densed matter research and its behaviour in magnetic field has been intensely investigated. 
The quality of graphene devices has improved over the recent years and ballistic transport 
over distances of several microns have been demonstrated recently [3, 4]. Since graphene is a 
gap-less semiconductor, it offers the possibility of creating internal interfaces with opposite 
charge carrier polarity. These so-called p-n interfaces [5-7] are formed by local electrostatic 
gating. 
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FIG. 1. Evolution of electron states in increasing magnetic field and design of a 
graphene p-n junction, a, A two-terminal graphene device consisting of a hole (blue) and an 
electron (red) cavity is sketched. By applying a weak field, the electron trajectories in the p- and 
n-cavities bend, leading to dispersing Fabry-Perot resonances, b, The field is increased until the 
cyclotron orbit becomes comparable to the cavity size, where resonant scar states can occur, c. The 
field is further increased and transport is still described by quasiclassical cyclotron orbits. Snake 
states along the p-n interface form, d. Finally, quantum Hall edge states propagate in opposite 
direction in the p- and n-region at higher fields, e, 3D design of the measured device. The Si02 
substrate is colored in blue and the bottom-gates in gold. The contacts, supported by the lift-off 
resist LOR (green) are colored in gray, f, Scanning electron microscope image of a device similar 
to the measured sample. The graphene is colored in blue and the bottom-gates in gold. The scale 
bar corresponds to 1 /rm. 

If electrons that propagate via skipping orbits enconnter snch a p-n interface, they will 
tnrn into snake states. These states consist of alternating half circles with opposite chirality 
and they transport cnrrent along the interface. Similar snake states have first been realized in 
GaAs/AlGaAs 2DEGs by defining regions of alternating magnetic field direction [8]. These 
states share the condition of commensnrability (similar to the above described magnetic 
focnsing experiments) with p-n snake states bnt they do not propagate along a single and 
tnnable interface. Snake states in graphene p-n jnnctions were claimed to have been observed 
in disordered snbstrate snpported samples [9] bnt the experiment lacked of direct evidence 
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for snaking trajectories. 

In this article we investigate ballistic transport across a graphene p-n junction in different 
magnetic field regimes and identify magneto-conductance oscillations as a direct signature of 
snake states. These findings are supported by detailed tight-binding simulations that allow 
us to visualize the alternating cyclotron orbits of the snake states. 


RESULTS 

Evolution of electron trajectories in graphene p-n junction 

In figure la-d we illustrate schematically how trajectories under increasing perpendicular 
magnetic held evolve in such a device. Figure la describes the low-held situation where 
transport is still dominated by Fabry-Perot oscillations with slightly bent trajectories [10]. 
As the held is increased resonant scar states (hg. lb) may occur, as observed in semiconduc¬ 
tor quantum dots [11]. At higher helds (hg. Ic) snake states at the p-n interface govern the 
electronic properties. Finally (hg. Id) the system enters the quantum Hall regime [12, 13] 
where transport is dominated by edge states and Landau level mixing at the p-n interface 
can occur [14, 15]. Even though this article is focusing on snake states, we will discuss the 
phenomenology of the mentioned magnetic held regimes. By doing so we present an integral 
picture of graphene p-n physics in magnetic held. 


Device architecture 

Figure le shows the design of the measured device and hgure If a scanning electron 
microscope picture of a similar sample. In a suspended 2x2 pm graphene sheet a p-n 
junction is formed by applying diherent voltages on the left (Vjeft) and right (Vj-ight) bottom- 
gates, resulting in diherent charge carrier concentrations nieft and Unght- The fabrication 
follows partly ref. [16] which we combined with a wet transfer process allowing us to align 
the graphene with the bottom-gates. For details see ref. [17] and Methods. 
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FIG. 2. Device characterization in the Fabry-Perot and in the quantum Hall regime. 

a, Two-terminal conductance as a function of left and right gate voltage shows regular Fabry-Perot 
oscillations at zero magnetic field. The inset reveals the narrow Dirac dip along the pp-nn diagonal 
(blue dashed line) from which a mobility of ^ ~ 470 • 10^ cm^V“^s“^ is deduced, b, Cuts along the 
same diagonal at different magnetic field strengths exhibit the expected quantum Hall plateaus at 
2, 6,10 e^//i. The G = 2 jh plateau is already visible at 60 mT. c, The color plot as a function of 
Ideft and Fright at 200 mT shows quantum Hall plateaus in the unipolar region, d, The numerical 
derivative dG/dV (in arbitrary unit) is recorded as a function of gate voltage V and B displaying 
the dispersion of the Fabry-Perot oscillations. V is the magnitude of gate voltage in the situation 
of antisymmetric charge densities nieft = —bright (red dashed line in a. i.e. the np-pn diagonal). 
The white dashed curve indicates the line along which the cyclotron radius Rc is equal to the cavity 
length L = 0.8 /rm, the region Rc < L is darkened and will be discussed in the main text, e, The 
measured pattern is reproduced by a tight-binding quantum transport calculation based on the 
designed geometry of the measured device. The small inset shows a resonant electron trajectory at 
low magnetic field. Constructive interference occurs if the phase along this trajectory is an integer 
of 27r, leading to the numerical solution of the yellow lines (see Methods). 
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Measurements in the Fabry-Perot and quantum Hall regime 


In the following we characterize the measured device in the zero, low- and high-held 
regimes. Figure 2a shows a two-dimensional color map of the electrical conductance 
(^(Kight, h^ieft)- As soon as a p-n interface is formed, G is lowered drastically. Regular 
Fabry-Perot resonances in the left/right cavity are visible as oscillations perpendicular to 
the zero density line in the left/right cavity (horizontal/vertical white dashed line in £g. 
2a), indicating ballistic transport [3, 4]. The inset in hgure 2a is a slice along the pp-nn 
diagonal (blue dashed). We estimate the mobility to be /x = ~ 470 ■ 10^ cm^V“^s“^ at 

carrier density n = 1.1 ■ 10® cm“^, which we calculated by a simple parallel plate capacitor 
model. The mobility is mainly limited by scattering at the contacts [3, 18]. 

In hgure 2b we show cuts along the pp-nn diagonal at diherent magnetic helds B and 
obtain quantum Hall plateaus aX G = Gq - where u = 2,6, 10,... is the hlling factor [12, 13] 
and Gq = ^ jh is the conductance quantum. Even at helds as low as 60 mT, the v = 2 
plateau is visible. The colorscale map in hgure 2c taken at 200 mT shows that plateaus 
develop in the unipolar region with conductance values given by the lowest number of edge 
modes in the left or right cavity, i.e. G = G'o'min(z/right, i^ieft) = 2, 6,10e®/h. This observation 
compares well to experiments of ref. [14], even though we apply only 0.2 T instead of 4 T. 
In the bipolar region the conductance stays well below 2 e®/h due to the smoothness of the 
p-n interface. 

In a further step we study the dispersion of the Fabry-Perot interference pattern in low 
magnetic held [10]. Figure 2d shows the numerical derivative dG/dV as a function of B 
and V, where V represents the magnitude of gate voltage in the situation of antisymmetric 
charge density (red dashed line in hg. 2a i.e. the np-pn diagonal). In this conhguration 
the device consists of two Fabry-Perot cavities of equal length L ^ 0.8 /xm. The darkened 
region where the cyclotron radius Rc < L will be discussed later. In hgure 2e a tight- 
binding transport calculation is shown which reproduces the measured interference pattern 
very well (for details see ref. [19] and Methods). We highlight the quality of the measured 
graphene and the ability of the simulation to capture the complex oscillation pattern of this 
micron-sized system in magnetic held. The dispersion of the Fabry-Perot oscillations can 
be described by bent electron trajectories such as the one sketched in the inset of hgure 2e 
[3, 10]. The condition for constructive interference is met if the accumulated phase along 
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such a trajectory is a multiple of 27r (see “Methods”). The yellow lines in hgure 2e are 
numerical solutions based on such a condition. 


Measurement of snake states 

We now discuss the regime where snake states emerge. In hgure 3a a snake state at a sharp 
p-n junction is sketched. Consider a charge carrier trajectory starting at the gray cross with 
momentum k in —x direction. Due to the magnetic held the trajectory is bent towards the 
p-n interface within the cyclotron diameter 2Rf. = 2hk/eB = 2fiy/nTT/eB. If the trajectory 
hits the p-n interface, the hole will be transmitted to the n side with high probability due 
to Klein tunneling [20]. At the upper edge of the sample, the snake trajectory scatters at 
the left side, resulting in a current towards the left contact. At lower n, Rc is reduced. As 
sketched in hgure 3b, the snake trajectory scatters to the right at the upper edge, resulting 
in a net current towards the right contact. With this mechanism one expects conductance 
oscillations that depend on B and n, and constant conductance along curves where Rc is 
constant. In hgure 3c we display calculated functions n{B) for constant Rc- Snake states 
occur once the cyclotron diameter 2Rc is smaller than the sample width W (green dashed 
curve) and can be described by quasiclassical trajectories as long as Rc is larger than the 
magnetic length Ib = \/h/eB (red dashed line). In the regime W > 2Rc > Ib additional 
parabolic lines show the condition for which the number of oscillations in the snake pattern 
is hxed and commensurate with W, i.e. m ■ 2Rc = W with m = 1, 2, 3,.... 

In hgure 3d we show the measured conductance G{B, V) (bottom) and its numerical 
derivative dG/dB (top). The measurement exhibits strong oscillations that follow parabola- 
like curves. We notice that the oscillations occur on a background of strongly decreasing 
conductance from G ~ 6e^/h to G < 2e^ jh. The steep decrease indicates that the transport 
becomes dominated by the low density region close to the p-n interface and this happens 
when 2Rc < W. 

In a real p-n interface the density does not sharply jump from the p to the n side but 
evolves smoothly. An electron trajectory in such a smooth p-n interface is sketched in hgure 
3e. Here the density gradient leads to an electric held E^, and its interplay with the perpen¬ 
dicular magnetic held results in the so-called ExB drift (here along y), leading to elongated 
cyclotron orbits. The condition Ry = const, can be studied in the measurement of hgure 
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FIG. 3. Parabolic-like conductance oscillations as a signature of snake states, a, Charge 
carrier trajectory (white) along a sharp p-n junction in perpendicular magnetic field starting at 
the gray cross, where Rc is the cyclotron radius and n{x) is the electron density, b, At lower p- 
and n-density Rc is reduced. In contrast to a, the trajectory results in current flow towards the 
right contact, c, Curves of constant Rc as a function of n and B. The continuous lines are given by 
the condition that the cyclotron diameter 2Rc is commensurate with the sample width W = 2 /rm. 
Snake states occur between the green dashed 2Rc = W and the red dashed Rc = Ib line. The 
black dashed line indicates up to which field transport is dominated by bent Fabry-Perot patterns 
(dark gray area). In the red area scar states can occur. Below the red dashed line Rc is smaller 
than the magnetic length Ib and Landau-levels start to dominate the transport (light gray area), 
d, Conductance as a function of antisymmetric gate tuning V and magnetic field is shown in the 
lower panel, and its derivative dG/dB (in arbitrary unit) in the upper panel. Striking lines of high 
and low conductance with a parabolic like B-dependence can be observed, e. In a smooth p-n 
junction (here: linear n{x)) the cyclotron orbits become elongated along the y-direction due to the 
additional electric field caused by the density gradient, f, dG(Fright, Fieft)/dF at 120 mT. The blue 
dashed lines in d and f are equivalent. The inset shows lines of constant electric field E^, at the 
p-n interface as a function of Fight, Feft taken from electrostatic simulations. 
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3f where we show (iG(Kight, at B = 120 mT. The measured interference pattern 

follows curves of constant E^. at the p-n interface (obtained from electrostatic simulations) 
as shown in the inset. 
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FIG. 4. Tight-binding transport calculation reproducing the experimental results and 
local current density profiles revealing the snake states, a, Tight-binding calculation of 
conductance G{B,V) (bottom) and its numerical derivative dG/dB (top, in arbitrary unit). The 
parabola-like lines seen in experiment are well captured, b, G(y) along the white dashed line in a 
at B = 90 mT. c, Calculated x-component of the local current density distribution, jx, for electrons 
injected from the left contact with a small DC-offset at i? = 90 mT and V = —4.8 V (dashed circle 
in b). The complex resonance pattern in the left cavity consists of “bubbling” trajectories that do 
not contribute to conductance. At the p-n interface (dashed line) an alternating current pointing 
to left (red) and right (blue) is observed, d-g. Current density profiles at different V corresponding 
to the circles in b. As a guide to the eye, a snaking trajectory following jx is added. From one 
conductance maximum to the next (i.e. c to g) one snake period is added. The snake state in e 
and f corresponds to a conductance minimum and c,d and g to a maximum. 
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Tight-binding theory of snake states 


So far we have seen that the oscillation pattern occurs in the regime where snake states 
are expected (i.e. 2Rc < W) and that the oscillations are related to transport along the 
p-n interface. We now present a quantitative comparison between experiment and theory. 
Figure 4a shows a quantum transport simulation of G{B,V) and dG{B,V)/dV based on 
a scalable tight-binding model [19] that fully takes into account the device geometry. The 
simulation compares very well to the measurement shown in hgure 3d. The parabola-like 
patterns are reproduced and a similarly steep decrease of conductance is obtained. In hgure 
4b a slice following the white dashed line in hgure 4a is shown. The visibility AG/G of 
the oscillations reaches 30% in theory and experiment and is enabled by the strong Klein 
collimation at the smooth p-n-interface. 

Next we apply the Keldysh-Green’s function method to extract local current density 
prohles (see Methods) at high and low conductance along this line. In hgure 4c we show the 
x-component of the current density jx, taken at V = —4.8 V (dashed circle in hg. 4b). The 
current is injected from the left contact using a small DC ohset. In the left cavity a complex 
resonance pattern appears, given by so-called ’’bubbling” trajectories [21] which are rehected 
before reaching the p-n interface and do not contribute to current between the contacts. The 
pattern relevant for transport is located at the p-n interface (dashed line). We observe that jx 
changes sign along the p-n-interface and that the blue and red regions penetrate the dashed 
line, indicating that transport is dominated by Klein-collimated snake trajectories. As a 
guide to the eye we added a curve in hgures 4d-g that follows the snake state. This is done 
for different V values for which G is maximal/minimal at a hxed magnetic held of 90 mT 
(colored circles in hg. 4b). In hgure 4d for example, the current density prohle corresponds 
to a conductance maximum where the current points to the right at the upper edge of the 
p-n interface. One period is added by changing V from one maximum to the next. More 
current density prohles evolving with V at hxed B = 90 mT are shown in Supplementary 
Movie 1. By tracing along one of the parabolic-like patterns, the current density prohle of 
the snake state stays constant (an example is given in Supplementary Movie 2). 

The conductance oscillates as a function of the ratio W/2Ry, the exact snake period 
corresponding to ARy is, however, difficult to determine using quasiclassical trajectories, 
since current is injected from many points under various angles resulting in a complex cusp 
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structure similar to what was predicted in refs [22-25]. The excellent agreement between 
measurement and calculated conductance for which we could determine local current density 
prohles clearly indicates that the oscillations result from snake state trajectories. 

There are additional parabola-like structures at lower magnetic field indicated by arrows 
in hgure 4a, these structures are however less pronounced in the experimental data. Those 
resonances occur in a regime where scar states (fig. lb) would be expected. The parabola- 
like behaviour indicates that the states are commensurate with a cavity dimension. In the 
model, the resonances disappear for non-reflective contacts as expected for scar states. 


DISCUSSION 

We investigated magneto-conductance of a ballistic graphene p-n-junction in the different 
magnetic held regimes. We have observed resonance patterns occurring in the intermediate 
quasi-classical regime in experiment and theory which result from the formation of snake 
states at the p-n interface. Among many other possibilities these states can be used to guide 
electrons on arbitrary paths with a high efficiency even at very low magnetic helds. This 
could be used to guide electrons away from sample edges to suppress uncontrolled momentum 
or spin-scattering. The directional scattering at the sample boundaries could be used to 
implement multi-terminal switches [24, 26]. Furthermore, the similarity between Andreev 
rehection and Klein tunneling is stressed in theory [27] leading to a correspondence of snake 
states and Andreev edge states which are of theoretical [28] and experimental [29] interest. 
Our work points out that snake states are highly tunable, occur at low helds and that 
ballistic graphene p-n-junctions in magnetic held reveal novel and intriguing phenomena. 
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METHODS 


Experimental methods. High quality graphene is obtained by in situ current annealing 

[30] . All measurements were performed in a variable temperature Helium cryostat with a 
base temperature of 1.5 K. We measured differential conductance G = dl/dV by standard 
lock-in technique applying an AC voltage of 0.1 mV at 77 Hz. 

For the quantum Hall data of hgure 2b-c we subtracted a contact resistance of 1.2 kO. 

We extracted the cavity length L used in hgure 2d-e from the spacing An between reso¬ 
nant Fabry-Perot peaks in the bipolar situation An = 2y^Fn/L and obtained L ^ 0.8 pm. 

Simulation methods. Real-space Green’s function method in the tight-binding framework 
using a scaled graphene Hamiltonian [19] is applied to simulate ballistic quantum transport in 
the present device, taking into account the realistic on-site energy prohle obtained by three- 
dimensional electrostatic simulation for the self-partial capacitances of the bottom gates. All 
the presented conductance simulations are obtained by calculating the transmission function 
at zero temperature, combined with the contact resistance 1.2 kfl. Local current densities 
are imaged by applying the Keldysh-Green’s function method in the linear response regime 

[31] based on the same model Hamiltonian used for conductance simulation. At each lattice 

site u, the bond charge current density J„ = is computed, where the 

summation runs over all the sites m nearest to n, is the unit vector pointing from n 

to m, and {Jn^m) is the quantum statistical average of the bond charge current operator 
Jn^m [32]. After computing for each site, the position-dependent current density prohle 
J{x,y) = [jx{x,y), jy{x,y)] is imaged. In £g. 4c-g the x-component jx{x,y) is shown. 

The low-held Fabry-Perot interference contours sketched in Fig. 2e are numerically ob¬ 
tained from solving the resonance condition A$ = 2j7r, arising from the path diherence be¬ 
tween the directly transmitted and twice rehected trajectories within the p cavity as sketched 
in the inset of Fig. 2e, which is found to be the major interference contribution. For such a 
simplihed model the phase diherence is given by A<F = <Fwkb + ‘^’ab + ‘^’Berry + *^*05 where 
*^’wKB = — {eBL/2hy-2L is the kinetic WKB-phase, <Fab = SA-eB/h is the Aharanov- 

Bohm phase due to the hux enclosed by the bent orbit segments, <FBerry = — vr(l — 
is the Berry phase, and $0 = tt is a constant phase due to rehections oh the two p-n junc¬ 
tion interfaces of the p cavity (smooth at middle and sharper at the contact side). Here 
kp is the numerical average of '\/7r[u(x]^'^7)y[ within the p cavity. The cavity size L is 
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numerically determined and is about half of the flake length Ll2 = 840 nm. The loop area 
is given by 6A = — sin0) with cj) = 2arcsin(L/2i?c) and Rc = hkp/eB. The form of 

the Berry phase follows from the consideration of [33] with the critical held estimated by 
Be = {2hkF/eL)\/l — Tc, where the critical transmission value is a parameter close to 
one and does not signihcantly inhuence the shape of the contours; Tc = 0.95 is chosen. The 
contours sketched in Fig. 2e correspond to j = 1,2, • • • ,8. 
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